Two-Center Integrals for r™- Polynomial Correlated Wave Functions 
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All integrals needed to evaluate the wave function of the form 



tot 



A 



L J 



EE 



for n = 1 and the Hamiltonian given are combined herein. For this 
form of the wave function, the integrals needed can be expressed as a 
product of integrals involving at most four electrons. An indication of 
how to increase or decrease exponents of r™- in steps of one or two is 
given. Some indication of how to proceed if the Hamiltonian contains 
1/rf - terms or if the wave function is of the form 



tot 



A 



Yl ">'>//, 

j <l < k 



is given 

r?Ai r mn with a < 



Consideration of all possible types of integrals involving 
b > 0, c > ; \a\ = \b\ = \c\ = 1, is given. 
Integrals are given in analytical form. These can be evaluated by nu- 
merical integration routines. 

INTRODUCTION 

Much success in ab initio calculations has been achieved with the use of 
correlated wave functions (wave functions that included the distance 
between two electrons explicitly). For the He atom 1 , the Li atom 2 , 
and the H2 molecule 3 , these wave functions have yielded the most ac- 
curate energy levels and molecular properties. Constructing the total 
wave function as a Slater determinant or antisymmetrized product is 
tantamount to using the Pauli principle, which excludes two electrons 
with identical quantum numbers and spin from occupying the same 
volume element at the same time. It does not tell us anything about 
two electrons with opposite spin, which we would expect to repel each 
other electrostatically. By including terms dependent upon the inter- 
electronic separation, we cause the probability, calculated from this 
wave function, of finding two electrons at specified regions of space to 
decrease when the two electrons approach one another. A correlated 
wave function can be an eigenfunction of spin and angular momentum. 
If .4(F) is an eigenfunction of the total and z component of spin and 
angular momentum, then ^ is an eigenfunction of the same 4 . Using 

1 



2 



the Nth-order permutation group 5 , all spin states 6 for the N electron 
system can be included: 

* = A{FR) 

F = F(l, 2, 3, . . . , TV) — (iV!)- 1 / 2 $ 1 (l)$ 2 (2) • • • 0> N (N) 

R = 1 + Yl + x ii r i + vn r % + • • • ) 

3 < I 

and A is the antisymmetrization operator. H is the Hamiltonian (in 
the Born-Oppenheimer approximation) in atomic units, Z a and Z& are 
the nuclear charges, R is the distance between nuclei and b, r iX is the 
distance between electron i and nucleus A , is the distance between 
electron i and electron j 

The coordinate system 7 used is confocal elliptical, is the out-of-plane 
angle and R is the distance between nuclei a and b. dri is the volume 
element and r i2 the interelectronic distance. We also have 

& = (r a i + r bi ) I R, Vi = ( r ai -r bi )/R 

1 < & < oo, -1 < Tji < 1, 

rfr, = i/? 3 (4 2 - d&drjidfa, < & < 2tt 

o 

^2 = ^ 2 {£l 2 + ^2 + Vi + % 2 - 2 - 2^x^77x772 

- 2 [(d 2 - l)(£f - 1)(1 - ^)(1 - % 2 )] V2 cos(0x - 2 )} (1) 

The basis functions $ s (j) and the total wave function are represented 
in Eqs. (2) and (3): 



= c;~i/^(C] - i) 7s/2 (i - v ]y°/ 2 e - a ^e^e im ^, 

s = (p 8 , Qs, Is, v s , a s , fi s , m,),i = V^l (2) 
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tot 



A 



L 3 



EE 

j < i 



(3) 



For molecules with more than two nuclei, the spherical coordinate sys- 
tem and Gaussian transforms 8 or ^-function expansions 9 can be used 
for integral evaluation. 

CLASSIFICATION OF INTEGRALS 

The classification of types of integrals involving Tl^r h kl r c mn can be consid- 
ered in the notation of picture writing 10 graph theory 11 . For the form 
of the wave function given in (3), the integrals needed can be expressed 
as a product of primitive integrals involving at most four electrons. All 
the integrals needed to evaluate this wave function involving to the 
first power are given. 



TWO - ELECTRON INTEGRALS 



We have 



(r 2 2 > = ($ s (2)r 2 2 $,(l)> = J rfr$ s (2)r 2 2 $,(l) 

i /»oo roc /*1 /*1 

x$ s (2)$i(l){[e i 2 + e2+^i+% 2 -2-266^2]5K;0)5(m t ;0) 

- [ (£i - lM - 1)(1 - Vl)(l ~ Vl)] V2 [ *K - 1; 0)5(m t + l; 0) 

+ S(m s + l;0)S(m t - 1;0)]} (4) 

The Neumann 12 expansion for 1 / r 12 in prolate elliptical coordinates is: 
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1=0 m= 



,21 + 1 



x Q\ m \&>i)Pl m \vi)P! ml (V2)e im ^-^ (5) 



(I- \m\)\ 
(/+ \m\)\ 

,|m| 



^ H (6< 2 ) 



\m\ . 



-P™(£) and are associated Legendre polynomials of the first and 

second kind in the conmplex plane 13 . £i<2 and ^2>i mean the smaller 
and the larger, respectively, of £i amd £ 2 - The -P™(£) and Q™(£) have 
the range [ 1 ,oo ]. The -P z m (^) have the range [ -1 , 1 ]. For further 



details, see Refs. 13 - 23. We have 

(& _ 1 W 2 Jl+m 



PT(v) 



2 l l\ d£ l+m 
{\-r) 2 ) m l 2 d l+1 



2 l l\ dr] l+rn 

QT(0 = (e 2 -i) m/2 ^Qi(0, 

<(i+l)/2 



Of 2 -!)'. 



- E 



2/ - 4j + 3 



Pi-2 j+ i(0 (6) 



The (I/V12) ,Eq. (10) can be expressed more concisely, using the defi- 
nitions of Eqs. (6)-(9), 

a = (p a , q a , 7a, v a , "a, fi a , m a ), b = (p 6 , g b , 7 b , i/ 6 , a 6 , /3 b , m b ), 

k,'A z ) = [ f d $Me - v 2 )e a v qa (e - i) W2 

x (1 - rfY a ' 2 e- a ^e + ^P^{i)P^{r]) (7) 



F°(z) = ^ 



[P°(z)]*(z* - 1) 



(8) 



z; = (-r 



(n - 









(9) 



(l/r 12 ) = ($ a (l)(l/r 12 )$ 6 (2)> = ivr^^ + m^O) 

x £ (2^+1)^ i^(*)^(*)*^. 6 (*)<fe (10) 



/i=(T=|m a | 



For I (7 a + a) and |(i/ a + a) integers, the one-dimensional integral 
K° ,^ a (z) can be evaluated analytically for each z and inserted in the 



numerical integration at the appropriate mesh points. 
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[{(i-a)/2] (-y a +a)/2[(n'-a)/2} (u a +a)/2 

k°^M= E E E E 

j=0 r=0 k=0 t=0 

(2//-2j)!(2//-2fc)! 
2»j\{ l i-j)\{ l x-a-2j)\2»'k\ 

(-l) [K+CT)/2] - r [|(7a + a)]\ (-1)«[§(i/« + a)]! 



X 



- k)\ {[i! — a — 2k)\ r\ [\{ la + a)-r]\t\ [\{v a + a)-t]\ 
X{4!Vo! ^^ < +1 (^2- S )! 

==0 u=0 a v 7 



5o!^EE 



s=0 u=0 



< +1 (S -s)\ 



h 



S = n - a + Pa - 2j + 2r, S 2 = S + 2, 

V = fi' -a + q a -2k + 2t, V 2 = V + 2 (11) 



The upper limit of j is §(// — cr) or — o — 1) , whichever is an integral. 
The upper limit of k is |(// — a) or |(// — <r — 1) , whichever is an 
integral. The upper limit of t is \{y a + cr) and the upper limit of r is 
|(7a+°") ; if these are not integrals, the summations are infinite ones. In 
practice, the , a (z) are evaluated recursively and numerically 24 ' 25 . 
If the $ s (j) of Eq. (10) are Slater-type orbitals, the integrals can be 
reexpressed 28 as a sum of "charge distributions". The r 12 expansion 
is needed for the evaluation of r\ 2 and for raising the value of n in r™ 2 . 
This expansion [ Eq. (12) ] has been derived by Harris 2 . The partial 
integration of Eq. (14) is used in the evaluation of the corresponding 
integral [Eq. (15)]. In Eq. (14), for z — > 00 and /i^O, the first term 
approaches 0. For yt, — , X° = 0; O per (fy denotes the interchange of 



a and b charge distributions. We have 



„ _ ^+2(6) | ^ 2 (6) , ig 2 fa) , P^M 

M if'(6) if'(6) p^m p^y 

" 4%) 4%) 4 CT| M 

ycr = (^-k| + l)(^-k|+2) 

" M (2^ + 3) 2 

" " (2/x-l) 2 

= v , 2(2/i + l)(4a 2 -l) 
" " (2/2- l) 2 (2/i + 3) 2 ' 

(/x-H)!2(2/i+l) 
" ~ (/i+H)!(2/i-l)(2/i + 3)' 

d^t^WP 1 [P^)] 2 (^ 2 -l)' 

= (2 M + % 2 - 2(/x - a + l)z^0^ - 1, (13) 



(r 12 ) = ($ a (l)$ 6 (2)r 12 > = —n 2 R 7 6(m a + m b ;0) 



+ 






+ (15) 



«(*): 



+ y; [ Kl_ 2 ^ a {z) + K^_ 2>a (z) ] + W;K^ a {z) (16) 
THREE-ELECTRON INTEGRALS 

The three - and four - electron integrals have been formulated in a 
straightforward manner, using partial integration. Equation (17) il- 
lustrates a technique of partial integration useful in the derivations. 
Whenever the product of two or more associated Lagendre polynomi- 
als occurs, these can be replaced by a sum over a single associated 
Legendre polynomial. The coefficients involve products of Clebsch - 
Gordon coefficients. Equation (22) is a Clebsch - Gordon series. For 
further information see Refs. 23 and 28 - 35. We have : 



f(y)dy g(t)dt 

h=y 



g(y)dy / f(t)dt, 



t=l 



u 



(l/) = f f(t)dt, v(y) = J\(s)ds (17) 



(ri 2 r 13 ) = ($ a (l)$ 6 (2)$ c (3)r 12 r 13 
1 



128 



rff'(5(m a + m b + m c ; 0) 



1 + 



per 



E E 

H=a=\mi,\ ji'=a'=\m c 



*{J dziSfy^W+XZSfy^z)] 

x + f;{z)n;-; a {z) k;^(z) 

+ X a »G°(z)N;^ a (z) K^ b {z)\ 

/OO ^ cr, cr' 

dz [F; (z)N^Jz) K^ b {z) 

+ F^(z)N^;l(z) K^ b (z) + X;g;(z)N$ : ;1(z) K^ b (z)} } 

(18) 

POO 

KM Z )= / F °M K l»M dz ( 19 ) 

J z 
poo 

KM Z )= / G;(z) K^ b (z) dz (20) 

J z 

POO 

KM Z )= / F;{z)K^ b {z)dz (21) 

J z 



i V j 

m m! (m + m') 



p\ m+m '\z) (22) 
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32 
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32 


3i 3 


mi 


m 2 


m 




m 2 


mi m 



X 



5(m; mi + m 2 ) 

(j - m)\ (ji + mi)! (j 2 + m 2 y 



(j + m)\ (ji - mi)! (j 2 - m 2 )\ 
x C(ji, j 2 ,j; 0,0,0) 

= 5(m;mi+m 2 )(2j + l) 1 ' 2 A(j x , j 2 , j) 
(-If 



C(ji,j 2 J;mi,m 2 ,m) 



x 



+ 32 - j - p)Kh -mi-p)\ 



[(j + m)\(j - m)\(ji + mi)!(ji - mi)!(j 2 + m 2 )\(j 2 - m 2 y] 1 / 2 
(j - j 2 + mi + p)\(j 2 + m 2 - p)\(j - ji -m 2 +p)\ 

A , , , \(a + b-cy.(b + c-a)\(c + a-b)\ ll/2 
A(a, b,c) = 



(a + b + c+ 1)! 

<£.(*) = /V - ^ 2 )^p;(0^'(0 

x P^(v)P^(v)^ Pa V qa (^ - l) 7a/2 (l - 7? 2 )^/ 2 e - a ^e +/3ar ' 
J J' - 

x{K^, a (z) + K^ a (z)} (23) 



/i 


\* 


f 




n 




J'~ 




a' 


m 




a 




m 



+ V a 
~ V 



ji — 2 jj,' J 



a m 



} 



j J' 
ii \i' J' 
a a' m 



Li + 2 fjf J 
a a' m 



{K?,j, a (z) + K™^ a (z)} 



= ^.Uc*) = ^ a + o E E { 



j j' 



( W 



- ( 



x ( u;, 



/i + 2 /i' + 2 J 

a a' m 

H + 2 /i'-2 J 

a a 1 m 

/i + 2 // J 

a a' m 

H // + 2 J 



1 v fx v fj/ 



H-2 n'-2 J 



rn 



li -2 n' + 2 J 



+ V a 



/i 



m 



a 

/i jJ, 



a 

2 // J 
i 



rn 



) 


/i 


V! 


J'~ 




a 


o' 


m 



a m 

' 2 J 

a a" m 

(*,*') 



} 



+ - K wa J N ZJ,a( z ) ( 25 ) 



(r 12 /r 13 ) = ($ a (l)$ 6 (2)$ c (3)r 12 /r 13 ) = — flV 
x5{m a +m b +m c ;0) ^ ^ (2//+l)Z;,'{ / dzF^z) 

/i=<r=|mf,| /i =<t =|m c | 



X ,c( z 



dzN; : ;, a (z) 



K;,,, b (z)^;^ c (z)G;(z) 

+K,^z)K°;^ jC (z)F°:(z)}} (26) 



10 



(ri 2 ri 3 /r 2 3) = ($ o(l) M 2 ) '* 'c(3)ri 2 r 13 /r ' 23 > 



R 10 n 3 5(m a + m b + m c ; 0) 
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1 + 0. 



per 



00/1 00 00 



^"=0 w"=—p," pJ=u' jl=(J 

x 5(iu; m 6 + «;")<?( M; cr) «5( | w"| ; 0(2/i" + l)Z;: 

/oo r ~ 



+ 



1: 



dz 



+ N^, a m^l(z)F^z) ] 



/oo — (<7,<r) 

/oo ^ 



X 



5/" 



-<t$y>( <7 > <7 ') 



: / (27) 



/■OO 

^«(*) =<£,.(*)= / ^{x)N$Jx)dx, (28) 
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P°Jx)P°',{x) 



X 



[P°{x)Y P$tx) 



-d 

dx 



x 



x 



(29) 



= %(*)=/ ^;,(x)Ar; : ;, a (x)rfx, (so) 



or: 



» 



/■oo 

/ <^^;;;>)^, (si) 



^ (32) 



: (<V) 
(/*>/*') 



a( Z ) = %a,uf)A Z ) 



E%$(x)Nffi^(x)dx, (33) 



^/i,//',a( Z ) 



=/ ^WV'.W^- (34) 



FOUR - ELECTRON INTEGRALS 
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We have 



(p*Dl) = ($ (l)$ 6 (2)$ c (3)$„(4) 
^R 13 7c 4 5(m a +m b +m c +m d ] 0) 



^12^13 ■ 

ri4 ' 



1 + 



per 



E E 

^=(7=1™;,! /i'=(7' = |m c 



x £ (2/ + 1)Z, 

^"=cr" = |m d | 



/oo 



X 



+ 



/oo 



x 



+ F^{z)K^^ 4 {z)^ b {z)] } 

(35) 



x 



L L' J J' 
J' 



/i /X 



J 



/X // J' 

<r <r' (<r + cr') 



(<t + <t') a" m 



Li" V 
(a + a') a" m 



a a' (a + a') 

{K?.L>a( Z )+ K ?>,L,a( Z )} 



(36) 



13 



= = s <'"• *•'•♦•*> E E E E 

L V J J 1 

'n n' J' 
a a' (a + a') 



X 



pi + 2 pi' J 
a a' (o- + o-') 

J pi" L] ' .V 

a + a') a" m 



pi -2 pi' J 
a a' (a + a') 

pi" L' 



} 



(7 + a') cr" m 

x + ^ >L »} + KKtvU z ) ( 37 ) 



~ (ct,<t',<t") ~ (<t',<t,<t") 

m ( ^>(*) = M <*w)» = E E E E 

J J' L U 



x{(u;u;; 



X 



/i + 2 /i' + 2 J 
a a' (a + a') 

pi + 2 pi' -2 J 
a 

J 1 



+V a V a , 



pi -2 pi' -2 



J 



a a' (a + a') 
pi- 2 pi' + 2 J 



PL pL 

a a' (a + a') 



a' (a + a') 

PL + 2 p! J 



a a' (tr + a") 
pL — 2 pi' 



a a' (a + a') 



+ V u 



a 



J 

a + a') 



x ( u;, 



'pi pi' + 2 J 



a a 



x 



J pi" L 
{a + a') a" m 



(a + a') 
J' 



pi pi 



' 2 



a 



J 

W + a') 



pi" V 



(a + a') p," m 



{KZ L , a {z)+K™, L (z)} 



(38) 



(r 23 ri 4 /ri 2 ) = ($ a (l)$ 6(2)$ c (3)$,(4)r 23 r 14 /r 12 > = ^R™n 4 5(a'; \m c \) 

CO CO CO 

x5(a"; \m d \)S(a; \m b +m c \)5(m a +m b +m c +m d ; 0) E ( 2 / i + 1 ) 

/i=cr \x'=c' fi"=cr" 

/CO ^ 
dzF;{z) [ ^W^W^+JV^^)^^) 

+ <£ 6 (*)*^,c(*) + ^sgftu*) + *KU*) ] 

(39) 



14 



i: 



dxF°„ (x) 



+K;:^,, d (x)N;: : ;jx)} (40) 



= I dxG;(x)K^ d (x)N; : ^ a (x) (41) 



(r 23 ri 2 /r 14 ) = ($ a (l)$ b (2)$ c (3)$,(4)r 23 r 12 /r 14 > = —R 13 n 4 5(a'; \m c \ 

olz 

CO OO CO 



X 



5(a"; \m d \)5(a; \m b +m c \)5(m a +m b +m c +m d ; 0) ^ ^ (V+ 1 ) 



[i=<j ii'=a' fj,"=cr" 



x{<;, b (^)^y, c (^)+iv ( ( ;;^ 

/•CO 



+ ^^), c , b W + ^(^')^W]) ( 42 ) 

5KLa(*) = /* w^(*)<$>)<k> («) 



VW*) = / W^to^&O*)^ (44) 



*W),c6(*) = / ^ 



(s) 



,c ( X ) ^M,/"' .6 ( X ) 

+^>)#S&(s)l **> (45) 



3KUaW=/ rf:(x)KZw(x)N*;^(x)dx, (46) 



KINETIC ENERGY INTEGRAL 
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Essentially no new basic integrals are involved in the evaluation of the 
kinetic energy, nuclear attraction, and overlap integrals. In some cases, 
a modified K a , s (z) integral is used. The modified integral H a , s (z) 
is defined in Eq. ( 47); it differs in that the £ 2 — rf term is not included. 



x e-^e^P^OPMd^, (47) 



-\j dT$ a (l)$ t (2)$ w (3)r{'sV? [r[ 2 <S> a (l)<S> x (2)<S> y (3)} = 
-ij(J + l) J rfr$ e (l)$ / (2)$,(3)r|' 3 rl2 2 

~ * [ dr Va( ^ 2) ^ (1) $/(2)$ g (3)rf3^ 2 , (48) 

J SI 71 
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$ a (l)$a(l) = $ e (l), ^(2)^(2) = $,(2), ^(3)^(3) = $ fl (3), 

D a (l) =p 2 a +p a + 2p a 'J a + la~ &1- ql- q a ~ 2q a Va ~ "a 

?2 



2a a & (Pa+la + 1) - 2fl,ifr (g a + i/ a + 1) + + ' ; " (/ " 



d 2 



^l 2 



2a aPa 2/3 a q a ^ 2 -gg fl2 ,,_i^[ "^-ffl 



K(l, 2) = (1 - ?? 2 )(g a + (3 aVl - - /3 a 6i»e0 

+ (?i - l)(Po - a ?i t h aa&mm) 

u 

- ^a^i + (f a - laj&mt.l'ni + 7af? 



+ «a£l + A,?7l +q a ~Pa , 2 , . 

x [fe 2 -l)(l-% 2 )te 2 -l)(l-^ 2 )] V2 



X 



J(4>l-<f>2) 



+ e 



-i(4>i—4>2) 



+m a (e 1 -vi) 



fe 2 -l)(l-^ 2 2 ) 
te 2 -l)(l-^i 2 ) 



1/2 



x |y(*i-*2) _ e -i(0i-^)j ; ( 49 ) 



For / = /' = , the kinetic energy integral [Eq. (48)] over electron 1 is 
the sum of — ^irRH^ e (oo) terms. For I = and I' — 1 , the integral is 
a sum of (ri 2 ) terms [Eq. (15)] ,evaluated using H°^, e (z) [Eq. (47)] for 
electron 1 and the usual K° , j(z) [Eq. (7)] for electron 2. For / = 
and /' = 2 , use [Eq. (4)] for (rf 2 ) with H a , {z) for electron 1. For 
/ = 1 and /' = , use the usual (l/ri 2 ) [Eq.'(lO)] and (r 12 ) [Eq. (15)] 
with [Eq. (47)] instead of [Eq. (7)] for electron 1. For / = /' = 1 , use 
(ri 3 /r 12 ) [Eq. (26)] and (r 12 r 13 ) [Eq. (18)] with [Eq. (47)] instead of 
[Eq. (7)] for electron 1. For I = 2 and V = , [Eq. (49)] is the sum of 



terms, 



1 

'2' 



-^7ri^ g (oo)(r? 2 > 



terms, and 
1 

~2"" v 4 
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terms. The (rf 2 ) are evaluated using [Eq. (47)] instead of [Eq.(7)] 
for electron 1. The case of the kinetic energy integral in which the 
Laplacian operates on r 12 and this is multiplied by r 12 is represented 
in [Eq. (50)]. 



1 



dr$ s (l)$ t (2)rf 2 V?K 2 $ a (l)$ :r (2)], (50) 



For I = V = 1 [Eq.(50)] equals 
A 



itR z Y8{m s + m a ; 0)5 (m t + m x ; 0)K$ Q e (oo)K° f (oo 



dr 



$/(2)r? 2 [D a (l) + £>.(l)]$ p (l) 



NUCLEAR - ELECTRON ATTRACTION AND OVERLAP 

INTEGRALS 

We have 

{Z a + Z b )Ci + (Z a - Z b )rn 



dr 



$ e (l)$ / (2)$,(3)ri 2 r 



i J.' 

13' 



(51) 



For I — I' — , the integral over electron 1 is 

- l -itR8{m e] 0)[(Z a + Z b )H% te (oo) + (Z a - Z b )H he (oo)] 

For V = and 1 = 1, use (r 12 ) [Eq. (15)] and [Eq. (47)] instead of [Eq. 
(7)] for electron 1. For I' — and 1 = 2, use (rf 2 ) [Eq. (4)] without 
the (^ 2 - r)l) term. For /' = / = 1 use (r 12 r 13 ) [Eq. (18)] with [Eq. 
(47)] instead of [Eq. (7)] for electron 1. The nuclear-nuclear repulsion 
integral is Z a Z b /R times the overlap integral [Eq. (52)] . 

(r[ 2 4 3 ) = [ ^(1)^(2)^(3)^3, (52) 



For I — I' — the integral over electron 1 is ^7iR 3 8(m e ; 0)Kq e (oo) . 
For I' — and / = 1 , the integral is (ri 2 ). [Eq. (15)] . For I' = and 
1 = 2, the integral is (rf 2 ) [Eq. (4)] . For V — I — 1 the overlap is 
(r 12 r 13 ) [Eq. (18)]. 

ADDITIONAL INTEGRALS 
Some integrals can be generated from previously given integrals by rais- 
ing or lowering the r^i index in even or odd steps. [Eqs. (53) and (54)] 
, using Eqns. (1) , (5) , (12), and (59). If the Hamiltonian contains the 
term 1/r?- , for the evaluation of spin-spin magnetic coupling or the rel- 
ativistic effects of an external electric field, then (1/rfj) [Eq. (55)] and 
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(1/rfj) [Eq. ( 56)] are some of the integrals needed. These results are 
based on a generalization of the Neumann expansion [Eg.(59)] 23 ' 31,36 ' 37 
. The C l n {x) [Eq. (58)] are Gegenbauer polynomials. For I = \ , the 
Gegenbauer polynomials are the same as Legendre polynomials. If the 
wave function [Eq. ( 3 )] is modified to be 



^tot = A 



1+ E EE^ 

j <l < k 



then terms (l/r^r^) [Eq. (60)] will occur in the kinetic energy inte- 
grals. 

(r 3 2 ) = J dr$ e (l)$/(2)T^ = \R 2 [ ($ (pe+2 ,... ) (l)r 12 $ / (2)) 

+ ($ (Pe ,, e+2 ,..)(l)r 12 $ / (2)) - 2($ e (l)r 12 $ / (2)> 

+ ($ e (l)r 12 $ (p/+2i ...)(2)) + ($ e (l)ri2$(p />ff/+ 2,..)(2)> 

- 2($ (pe+1 , ge+li ... ) (l)r 12 $ (p/+1 , g/+li ... ) (2)) - (l)r 12 $ / _(2)> 

-($ e _(l)r 12 * /+ (2)>], 

e+ = (Pe, Qe, 7e + 1, + 1, «e, m e + l), 

e+ = (p e , 5e, 7e + 1, + 1, «e, As, ™e ~ 1) (53) 



(rls/rn) = J rfr$ a (l)$ s (2)$ t (3)r 2 3 /r 12 = ^ { <<t> t (3)) 

x [ ($ (Pa+2i ... ) (l)(l/r 12 )$ s (2)) + ($ (Pai9a+2 ...)(l)(l/r 12 )<l» s (2)) 
- ($ a (l)(l/r 12 )$.(2)> ] + ($ a (l)(l/r 12 )$.(2)> 

x [ ($( K+2 ,..)(3)> + (<f( Ptl9t+2 ,...)(3)) - (* t (3)> ] 
-($ t _(3))($ a+ (l)(l/r 12 )$ s (2))-($ t+ (3))($ a _(l)(l/r 12 )$ s (2)) 
-2($ (pt+1)9t+1 .... ) (3))($ (p 

a+l,?o+l,- 

)(l)(l/r 12 )$ s (2))}, (54) 



(l/r 3 2 ) = J rfr$ a (l)$ 6 (2)(l/r 3 2 ) = 7r 2 i? 3 5(m a + m b ; 0) 

oo 2 1 

Z=0 m=|m a |+l 

/oo 
Fr(z)K™ a (z)K™ b (z)dz, (55) 
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(l/r 2 12 ) = J dr$ a (l)$ 6 (2)(l/r 1 2 2 ) = ^ 2 R 4 5(m a + m b ; 0) 



oo n 



n=0 Z=|m a | 



i + (-iy+ m « 



(l!) 2 (w-Z)!(n+l)(2l + l) 
(ra + Z)! 



x 



(Z - m a )!(Z + m a )! Z" 00 ^!An-W0^ L ^V*,H0^ 



a(0 = (Pa, Qa, 7a + Z, i/ a + Z, a a , /3 a , m a ), 

6(Z) = (p 6 , g b , 7 fe + I, v b + /, a 6 , /3 6 , ra 6 ), (56) 



x e^e^C l n (OC l Av), (57) 



[n/2] 



2 ™-2j(_l)j(7 + n _ j _ 1)! x n-2j 

j!(n-2j)!(i-l)! 



(58) 



The upper limit of the sum over j is \n or |(n — 1) , whichever is 
integral. We have 36 ' 37 



2r 



12 



-2p oo n 



J2J2 d M K 1 - tfx 1 - % 2 )(£i - ^ - !)] 



1/2 



n=0 Z=0 



X 



^(6>2)crk6<i)^(^)crk^)cr i/2 [cos(0 1 - 



p>0, p^- 
-2 2 ' +1 r(2p- l)[r(p + /)] 2 (n-Z)!(n+p)(2p + 2Z- 1) 



^(0 = -c, 



U J0 f 



[r( P )] 2 r(2p + n + /) 
(x 2 - \y u ' l / 2 dx 



d_ 



.cm. 



(£2 _ !)-«-l/2 

raop ' (59) 
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(l/ri 2 r 13 ) = J rfr$ a (l)$ 6 (2)$ c (3)l/r 12 r 13 = 
^7r 3 R 7 5(m a + m b + m c ; 0)6 (a; \m b \)6(a'; \m c \) 

oo oo 

x^^(2l + l)(2j + l)^;' 

l=a j=a' 

/oo 
dzN^(z)Kr^(z)^z)Fr^), (60) 
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